Mon. Not. R. Astron. Soc. OOO,[TH5](2013) Printed 26 February 2013 (MN MT^X style file v2.2) 



Further Evidence for the Accretion Disk Origination of the 
Double-Peaked Broad Ua of 3C390.3 



o 
in 

(N 



o 

u 

o 



> 
oo 

o 

O 

\o 

o 

m 



Xue-Guang Zhang^'^ 

^Purple Mountain Observatory, Chinese Academy of Sciences, 2 Beijing XiLu, Nanjing, JiangSu, 210008, P. R. China 
'^Chinese Center for Antarctic Astronomy, Nanjing, JiangSu, 210008, P. R. China 



ABSTRACT 

In the letter, under the widely accepted theoretical accretion disk model for the double- 
peaked emitter 3C390.3, the extended disk-like BLR can be well split into ten rings, 
and then the time lags between the lines from the rings and the continuum emission 
are estimated, based on the observed spectra around 1995. We can find one much 
strong correlation between the determined time lags (in unit of light-day) and the flux 
weighted radii (in unit of Rq) of the rings, which is well consistent with the expected 
results through the theoretical accretion disk model. Moreover, through the strong 
correlation, the black hole masses of 3C390.3 are independently estimated as ~ IO^Mq, 
the same as the reported black hole masses in the literature. The consistencies provide 
further evidence to strongly support the accretion disk origination of the double- 
peaked broad balmer lines of 3C390.3. 
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1 INTRODUCTION 

As one well-known double-peaked emitter and one well- 
studied mapped AGN, 3C390.3 has been studied for more 
than four decades (Burbidge & Burbidge 1971, Chen & 
Halpern 1989, Dietrich et al. 1998, 2012, Eracleous et al. 
1995, Eracleous & Halpern 1994, 2003, Flohic & Eracleous 
2008, Gezari et al. 2007, Gliozzi et al. 2006, 2011, Kol- 
latschny & Zetzl 2011, Leighly et al. 1997, Lewis & Era- 
cleous 2006, Netzer 1982, O'Brien et al. 1998, Perez et al. 
1988, Popovic et al. 2011, Sambruna et al. 2009, Sergeev et 
al. 2002, 2011, Shapovalova et al. 2001, 2010, Zhang 2011a, 
2013, Zu et al. 2011). Based on the long-term variabilities 
of the double-peaked balmer lines and the continuum emis- 
sion, the proposed double-stream model (Vellieux & Zheng 
1991, Zheng 1996) has been ruled out for 3C390.3 due to 
the un-obscured emitting region on the receding jet (Livio 
& Xu 1997), and the binary black hole model (Begelman et 
al. 1980, Boroson & Lauer 2009, Gaskell 1996, Zhang et al. 
2007) has been ruled out for 3C390.3 due to the unreason- 
ably large central BH masses (Eracleous et al. 1997), the 
proposed accretion disk model (see references listed above) 
has been widely accepted for 3C390.3. 

Furthermore, besides the analysis for the double-peaked 
line profiles, based on the reverberation mapping technique 
(Blandford & Mckee 1982, Home et al. 2004, Peterson 1993), 
some simple information of the dominant gas motions in the 
BLR of 3C390.3 has been estimated through the response 
of different parts of the double-peaked broad emission lines. 



e.g., the line core, peaks and wings (Dietrich et al. 1998, 
2012, Popovic et al. 2011, Sergeev et al. 2002, Shapovalova 
et al. 2001, 2010): the blue and red peaks of the double- 
peaked broad balmer lines vary with the same time delay 
relative to the continuum variations, which is well consistent 
with the expected results under the accretion disk model 
for the double-peaked broad lines. Furthermore, the results 
in the literature have shown that there are no time delays 
between line wings, peaks and line core of double-peaked 
broad balmer lines of 3C390.3. 

However, besides the simple results above based on the 
reverberation mapping technique applied to the line wings, 
peaks and core, there is one another interesting way to check 
the proposed accretion disk model for 3C390.3. Through 
the widely accepted accretion disk model, the proposed ex- 
tended disk-like BLR of 3C390.3 can be well split into mul- 
tiple rings with different radii. Apparently, time lag between 
one line from one of the rings and the continuum emission 
should sensitively depends on the radius of the ring. More- 
over, the time lags are much different from the ones in the 
literature for the line wings, peaks and core, because the 
line core includes great contributions from the line photons 
coming from all the rings of the disk-like BLR of 3C390.3 
(such as the following results shown in Figure 1). Certainly, 
our results are much different from the results in Flohic & 
Eracleous (2008): each observed line profile of 3C390.3 are 
split into 36 parts, each part has one narrower wavelength 
range ~ 20A. However, in our letter, each observed line is 
split into ten parts, each part from the corresponding ring 
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of the BLR has the full wavelength range of the observed 
Ha. To check the dependence of the time lags on the radii 
of the rings is our main objective of the letter, which could 
provide further evidence to support or against the accretion 
disk model for the double-peaked emitter 3C390.3. 



2 MAIN RESULTS 

As shown in our previous papers (Zhang 2011a, 2013, Pa- 
per I and Paper II), we have shown that the standard el- 
liptical accretion disk model with few effects from prob- 
able bright spots and/or warped structures (Eracleous et 
al. 1995) can be well applied to describe the properties 
of the observed double-peaked broad Ha of 3C390.3 ob- 
served around 1995. So that, in the letter, the standard 
elliptical accretion disk model is mainly considered, and 
there are no further considerations for the other models 
which can be well applied to 3C390.3 observed in differ- 
ent periods. And moreover, because the broad Ha hav- 
ing stronger intensity and less contamination from other 
narrow emission lines, we mainly consider the 67 spectra 
with available broad Ha from the AGN WATCH project 
(http:/ /www. astronomy.ohio-state.edu/~agnwatch/ ) (Diet- 
rich et al. 1998). 

As we simply discussed in the Introduction, it is very 
interesting to check the time lags between the continuum 
emission and the lines from the rings with different radii, 
under the accretion disk model for 3C390.3. Here, the to- 
tal proposed disk-like BLR of 3C390.3 under the stan- 
dard elliptical accretion disk model is evenly split into ten 
rings (i = 1,2, 3. ..10) with lower and higher boundaries, 
[ro + S{r) X i,ro + 5{r) x (i -I- 1)], where ro and ri represent 
the inner and outer boundaries of the total BLR, 5{r) = 
(ri — ro)/10 means the extended size of each ring. Based 
on the best fitted results for the 67 observed double-peaked 
broad Ha, the mean value of ro is around ~ 200Rg and the 
mean value of ri is around ~ 12QQRg, where Rg = '^^^^^ 
is the Schwarzschild radius. In other words, each observed 
double-peaked broad Ho? could be well split into ten parts 
(ten lines), and each part has apparent double-peaked line 
profile. Here, the BLR would not be split into more rings, 
otherwise light traveling time for each ring should be much 
shorter than the average time interval for the light curves of 
3C390.3. After the consideration of the more recent black 
hole masses of 3C390.3 Mbh ~ lO^M© (Dietrich et al. 
2012), the light traveling time for each ring with extended 
size of ~ IOORg is about 6days, which is one appropriate 
time value. Figure 1 shows one example for the ten lines from 
the ten rings. It is clear that the line core defined in the lit- 
erature includes apparent and strong contributions from the 
lines coming from all the rings. Therefore, to check the time 
lags between the continuum emission and the lines from the 
rings should be more meaningful. Then the intensities for 
the ten lines from the ten rings, fluxi{i = 1, 2..., 10), can be 
well determined by the accretion disk model. Certainly, we 
should note the values of fluxi could not be directly used, 
before some procedures being applied to complete the fiux 
calibration. 

There are three steps being applied to complete the 
flux calibration for the line intensities of the ten lines from 
the ten rings. First and foremost, because the 67 spectra of 
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Figure 1. Example for the ten lines coming from the ten rings 
of the disk-like BLR of 3C390.3. Top panel shows the clear pro- 
files of the ten lines from the ten rings, difi'erent colors and line 
styles represent lines from diflFerent rings. Color of red, green, 
blue, yellow, and line style of solid line, dotted line and dashed 
line show the line from the ring with increasing radius. Bottom 
panel shows the sum of the ten lines (thick solid line in red color) , 
which further represents the best fitted results for the observed 
double-peaked Ha around JD- 2449870. Moreover, the shadow ar- 
eas in the bottom panel show the masked regions for the narrow 
lines. 

3C39Q.3 around 1995 are observed by different instruments 
in different observatories under different configurations, the 
measured line intensity based on the observed line profile 
should have low confidence level. Therefore, we collected 
the reliable values {Fagn with corresponding uncertainties) 
listed in the AGN WATCH project, with the completed fiux 
calibration having been applied (Dietrich et al. 1998) for 
the observed lines. Besides, the line intensities in the AGN- 
WATCH project include contributions from the narrow lines 
around Ha. In order to ignore the effects of narrow lines on 
the following determined time lags, the intensities of the 
narrow lines should be subtracted from Fagn- Here, based 
on the intensity of [OIII]A5QQ7A in Dietrich et al. (1998) 
and the fiux ratio of [Nil] doublet and narrow Ha to [QUI] 
in Dietrich et al. (2012), the total intensity of narrow lines 
around Ha is determined as 

Fnarrow = flux{Han) + flux{[NII]) 

^ (3.53 + 1.34 + 0.44) x flux{[OIII]) (1) 

~ 117 X 10"^^erg/s/cm^ 

Finally, the available line intensities [Fluxiii — 
1, 2, 3, .., 10)) of the ten lines from the ten rings for 
each observed broad Ha can be determined by the inten- 
sities (fluxi) from the direct observed lines without fiux 
calibration, 

Flux, = J}^^' X {Fagn - F„arrow) (2) 
2^ J lux i 

And, the uncertainty of FluXi can be determined by 
the uncertainties of Fagn and flux{[OIII]). Then, 
the light curves of the continuum and the ten lines 
from the ten rings can be determined and shown in 
Figure 2. Here, the light curve of the continuum is 
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Figure 2. The results based on the SPEAR method and the CCF method applied to the light curves of the continuum emission and 
the ten lines from the ten rings. For each line and continuum, there are four panels, two panels show the light curves of the line and the 
continuum (red dots with error bars) and the the best descriptions for the light curves by the SPEAR method (green lines), one panel 
shows the distributions of the expected time lags between the line from the ith ring of BLR and the continuum through the SPEAR 
method (the distributions is normalized to have a maximum value equal to the peak value), and one panel shows the CCF results for the 
two light curves. The bottom panel shows the corresponding results for the continuum emission and the total observed broad Ha. In the 
figure, the vertical dotted lines show the positions with zero time lags, the blue solid lines show the peak positions of the distributions 
of the time lags through the SPEAR method, and the peak positions of the CCF results. The digital numbers in blue color gives the 
corresponding time lags. In the panels for the light curves, the unit for the continuum is 10~^^erg/s/cm^/A, and the unit for the emission 
lines is 10~-^^erg/s/cm^ 



the well-time binned one from the AGNWATCH Project 
(http:/ /www.astronomy.ohio-state.edu/^agnwatch/ScSQO/lcv/ ). 

Now, based on the light curves, we can check the depen- 
dence of the time lags for the ten lines on the radii of the 
ten rings. There are several commonly accepted CCF (Cross 
Correlation Function) methods to estimate time lag between 



two data series, such as the ZDCF method (Z-transfer Dis- 
crete Correlation Function, Edelson & Krolik 1988, Peterson 
1993, White & Peterson 1994), the ICCF method (Interpo- 
lated Cross Correlation Function, Gaskell & Peterson 1987, 
Peterson 1993), the MCCF method (Modified Cross Cor- 
relation Function, Koratkar & Gaskell 1989, Koptelova et 
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al. 2006) etc.. More recently, Zu et al. (2011) have reported 
one new method to estimate time lag between continuum 
and emission line of AGN (SPEAR method, Stochastic Pro- 
cess Estimation for AGN Reverberation), based on the as- 
sumption that all emission-line light curves are time-delayed, 
scaled, smoothed, and displaced versions of the continuum. 
This alternative approach fits the light curves directly us- 
ing a damped random walk model (DRW model, Kelly et 
al. 2009, Kozlowski et al. 2010, MacLeod et al. 2010) and 
aligns them to recover the time lag and its statistical confi- 
dence limits. In the letter, both the ICCF method and the 
Zu's SPEAR method are applied to estimate the time lags 
between the continuum emission and the ten lines from the 
rings. 

When the ICCF method is applied, in order to ignore 
the efi^ects of the large time gaps of the light curve as far as 
possible, only the data points within JD from JD-2449700 to 
JD-2450000 (around SOOdays) are considered. Then the cor- 
responding uncertainty of the time lag can be determined by 
the Monte Carlo method as what we have done in Zhang et 
al. (2011b). When the SPEAR method is applied, in order to 
obtain good enough results about the light curves through 
the DRW method, all the data points in the light curves are 
considered. Here, we do not describe the SPEAR method in 
detail, which can be found in Zu et al. (2011) and in the web- 
site: https://bitbucket.org/nyel7/javelin. When the SPEAR 
method is applied, the number of walkers, burn-in iterations 
and sampling iterations for each walker are 300, 300 and 
300 respectively in the MCMC (Markov chain Monte Carlo 
method) analysis. Moreover, when the SPEAR method is 
applied, range from -lOdays to lOOdays is set as the bound- 
aries for the time lag. The final results based on the ICCF 
method and the SPEAR method are shown in Figure 2. 

Based on the results shown in Figure 2 (especially the 
distributions of the time lags through the SPEAR method, 
and the CCF results), the time lags between the continuum 
emission and the ten lines from the ten rings are increasing 
with encreasing the radii of the rings. More clearer results 
are shown in Figure 3. In Figure 3, we show the correlations 
between the flux weighted radii of the ten rings (in unit of 
Rg) and the time lags determined by the SPEAR method 
and the CCF method. Here, the flux weighted radii of the 
ten rings with lower and higher boundaries are calculated 
through the standard elliptical accretion disk model. For 
the ith (i — 1,2, 3. ..10) ring, there are 67 flux weighted 
radii based on the lower and higher boundaries for the 67 
rings of the 67 observed broad Ha. Then, the mean value is 
accepted as the flux weighted radius of the ith ring, and the 
corresponding uncertainty is estimated as the value of the 
flux weight radius of the ith ring minus the minimum inner 
boundary for the 67 ith rings of the 67 observed spectra. It 
is clear that there are well consistent time lags determined 
by the SPEAR method and by the CCF method. Moreover, 
we can find the time lags between the ten lines and the 
continuum emission are strongly linearly correlated with the 
radii of the ten rings. 



^ = (19.45 ± 2.15) X , '^fP'^^^ 
Rg fight — days 



(16.87 ± 1.91 X 



TCCF 



(3) 





Figure 3. On the correlation between the time lags determined 
by the SPEAR method and by the CCF method (left panel), and 
the correlation between the time lags and the radii of the rings 
(right panel). In the left panel, the sold circle represents the values 
for the total observed Ha. In the right panel, circles are for the 
time lags based on the SPEAR method, triangles are for values 
based on the CCF method, solid circle and solid triangle are for 
the total observed Ha. In the left panel, the solid line represents 
TCCF = '''SPEAR- Iri the right panel, the solid line and the dotted 
line are the best fitted results (Equation 3) for the correlation 
about TspEAR S'Ud for the correlation about tqcf- 



unit of Rg- The corresponding values of Chi-square divided 
by number of degrees of freedom are 0.04 and 0.4 for results 
about TSPEAR and tccf- The results above are strongly 
support the accretion disk model for 3C390.3. 



3 DISCUSSIONS AND CONCLUSIONS 

Before further discussion, we firstly check whether our re- 
sults about the time lags are reliable, especially for the time 
lag between the total observed broad Ha and the contin- 
uum emission. In Zu et al. (2011), the time lag between the 
continuum emission and the H/3 determined by the SPEAR 
method is about 27.9l^ g light-days for 3C390.3, which is 
well consistent with our result of ~ 25 ± 5 light-days be- 
tween the continuum emission and the observed Ha by the 
SPEAR method. Moreover, the time lag between the contin- 
uum and the observed Ha is about 20 ± 9 light-days by the 
CCF method in Dietrich et al. (1998), which is consistent 
with our result 27 ± 5 light-days by the CCF method. The 
consistencies above indicate our procedures to estimate the 
time lags by the SPEAR method and by the CCF method 
are reliable. And therefore, the results shown in Figure 2 
and in Figure 3 have high confidence levels. In other words, 
the strong positive correlation shown in Figure 3 are reli- 
able, which strongly support the accretion disk model for 
3C390.3. 

Besides to provide strong evidence for the accretion 
disk model for 3C390.3, the strong positive linear correla- 
tion shown in Figure 3 could provide one another way to 
independently estimate the black hole masses of 3C390.3, 



light — days 

where R means the flux weighted radii of the ten rings in 



Mbh 



IOSMq 0.005687 x ^/jr 



(4) 
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where the factor of '0.005687' is the value used to trans- 
fer the unit of Rg to the physical unit of light-days. Then, 
based on the Equation 3, the black hole masses of 3C390.3 
are around 9.3l?;? x 10® Mq if R oc TSPEAR being considered, 
and the black hole masses are around lOAt\ l x lO^M© if 
R oc TccF being considered. The black hole masses are very 
well consistent with the more recent result in Dietrich et 
al. (2012): Mbh ~ lO'-^M©. The results not only indicate 
our results in Figure 3 are reliable, but also provide one op- 
tional independent method to determine black hole masses 
for double-peaked AGN besides the M-sigma method (Fer- 
rarese & Merritt 2001, Gebhardt et al. 2000, Giiltekin et al. 
2009, Lewis & Eracleous 2006) and the Virialization method 
(Dietrich et al. 2012, Peterson et al. 2004, Woo et al. 2010). 

Finally, we give our summary as follows. The time lags 
between the ten lines coming from the well split ten rings of 
the disk-like BLR and the continuum emission have been cal- 
culated through the SPEAR method and the CCF method, 
under the theoretical accretion disk model. Then, we can 
find one strong correlation between the time lags and the 
radii of the rings. Moreover, based on the correlation, the 
independently determined black hole masses of 3C390.3 are 
well consistent with the more recent values in the literature. 
The results above give further and strong evidence for the 
accretion disk model for 3C390.3. 
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